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AMENDED CLAIMS 



1. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
digital signal defined by a series of N real starting 
samples x(n)/ with N a power of two and n e [0..N-1]/ 
comprising successive transformation sj:eps (2) for 
transforming input samples into output samples,^ all the 
transformation steps being performed by means of a 
single set of butterfli es with several inputs and 



10 several outputs, the operating mode of which is ^ i^^^ 
modified selectively in each transformation step, the 
input and output samples of each transformation step 
being stored in a storage memory, a series of N output 
samples y(n) representative of the fast Fourier 

15 transform or the inverse fast Fourier transform of the 



yj 



output samples x(n) being provided in the last 
transformation step, 

characterized in that output samples y(n) are 

real, 

20 and in that the output samples of a butterfly 

replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples x(n) processed in the first transformation step 
are classified in bit-reversed order of their index n, 

25 output samples y(n) are provided in the last 
transformation step in ascending order of index n, 
these output samples being defined by the following 
relations : 

y (0) = Re[X(0) ] 



y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-1 

y(n) = Im[X(n/2)] for* n being even and 

different from 0 

y(N-l) - Re[X(N/2) ] 
where samples X(n), with n € [0..N-1], designate the 
complex samples of the series corresponding to the fast 
or inverse fast Fourier transform of the starting 
sample series x (n) . 

2. A method of calculating the fast Fourier 
transform or the inverse fast Fourier transform of a 
digital signal defined by a^^eries of N complex ^am ples 

X(n) conjugated _Jby pairs, characterized in that the 

calculation is jdone on a series of N real starTTng 

samples yjn) representative of the series o f com plex 

sainples X(n)^,_ with N power of two and n g [0.,N-1], 
the starting samples y(n) being defined as follows: 

y(0) = Re[X(0)] 

y(n) = Re [X ( (n+1) /2) ] for n being odd and 

different from N-1 
y(n) = Im[X(n/2)] for n being even and 

different from 0 

y(N-l) = Re[X(N/2) ] 
in that this method comprises successive transformation 
steps for transforming input samples into output 
samples, a series of N real output samples x(n) 
representative of this fast or inverse fast Fourier 
transform being provided in the last transformation 
step, all the transformation steps being performed by 
means of a single set of butterflies with several 
inputs and several outputs, the operating mode of which 
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is modified selectively in each transformation step, 
the input and output samples of each transformation 
step being stored in a storage memory, 

and in that the output samples of a butterfly 
replace the corresponding input samples of the same 
rank in the storage memory, so that, if the starting 
samples y(n) processed in the first transformation step 
are classified in ascending order of index n, the 
output samples x (n) are provided in the last 
transformation step in bit-reversed order of index n, 

3. The calculation method according to claim 1 or 
characterized in that, in each transformation step, 

each butterfly transforms input sample pairs, the ranks 
of the input samples of the same pair within the series 
of input samples of said transformation step being 
symmetrical with respect to a center between the end 
rank values of the input samples transformed by said 
butterfly, 

4. The calculation method according to claim 3, 
characterized in that it comprises ii-1 transformation 
steps Ep with u=log2(N) and p g [0..vi-2], 

5. The calculation method according to claim 4, |in 
turn dependent on claim 3, in turn dependent on claim 
l,j characterized in further comprising: 

a preliminary step of modifying the sequence of 
the starting samples x(n) ranked in ascending order of 
index n and showing them in bit-reversed order of index 
n in the first transformation step, and 

a final step of processing the series of output 
samples y(n) and providing a series of N complex 
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conjugated samples X(n) corresponding to the fast or 
the inverse fast Fourier transform of the series of 
starting samples x{n), 

6. The calculation method of claim 4, (in turn 
5 dependent on claim 3, in turn dependent on claim 1, or 
according to claim 5, J characterized in that, in each 
transformation step Ep, butterflies are distributed 
among N/2P^^ calculation blocks, 

in that each calculation block has a peripheral 

10 butterfly and/or 2^-1 internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 
the input samples of rank 2^^^a, 2P^^a+2P^'-l, l^^^a-^l^^^ 
2^^^a+2^"^^~l into output samples of the same rank, 

15 and in that the internal rank x butterfly of the 

rank a calculation block in transformation step Ep 
transforms the input samples of rank 2^''^a+2T+l, 
2P^W2T+2, 2P^^a+2P^^-2T-3, 2P^W2P^2-2x-2 into output 
samples of the same rank, with P>1 . 

20 7. The calculation method according to claim 6, 

characterized in that each butterfly is assigned a 
coefficient W^, whereon the calculation inside the 
butterfly is based, said coefficient being equal to e" 
j(27cs/N) ^ith s € [0.. N/4-1] for a fast Fourier transform 

25 and is equal to e^*^"^^''* with s g [0..N/4-1] for an 
inverse fast Fourier transform. 

8. Calculation method according to claim 7, 
characterized in that the internal rank x butterfly of 
the rank a calculation block in transformation step Ep 

30 is assigned coefficient with 5 - (x+1) . (N/2P'"^) . 
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9. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
have 

5 - four inputs for receiving input samples and 

four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
inputs, 

10 in order to selectively apply different 

transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
and a coefficient admitted on said corresponding 

15 additional inputs. 

10. The calculation method according to claim 9, 
characterized in that, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 
an internal butterfly, 

20 in that the permutation signal is 0 for the even 

rank calculation blocks, including rank 0, and 1 for 
the other ones. 

11. The calculation method according to claim 10, 
characterized in that, in transformation step Ep, each 

25 calculation block comprises one peripheral butterfly 
and 2^-1 internal butterflies. 

12. The calculation method according to claim 11, 
characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the last 

30 transformation step, and otherwise 0. 
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13, The calculation method according to claim 12, 
characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j.B, the 
butterfly delivers the following output samples si, s2, 

5 s3, and s4 

1) if the primary mode and secondary mode signals 
are 0: si = el + e2 

s2 = el - e2 
s3 = e4 - e3 
10 s4 = e3 + e4 

2) if the primary mode signal is 0 and the 
secondary mode signal is 1: 

si = el + e2 + e3 + e4 
s2 - el - e2 
15 s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B. e4 
20 s2 = e2 + B.e3 + A.e4 

s3 = el - A.e3 + B.e4 
s4 = -e2 + B,e3 + A,e4 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 

25 si = el - A.e3 + B.e4 

s2 = -e2 + B.e3 + A.e4 
s3 = el - A.e3 - B,e4 
s4 = e2 + B.e3 + A.e4 

14. The calculation method according to claim 10, 
30 characterized in that, in transformation step Ep,. each 

calculation block comprises: 
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2^-1 internal butterflies and a peripheral 
butterfly for the even values of index p as well as for 
the last transformation step if p is even, and 
2^-1 internal butterflies, otherwise, 
5 15, The calculation method according to claim 13, 

characterized in that the secondary mode signal is 1 if 
the peripheral butterfly is used for the last 
transformation step with p being odd, and otherwise 0, 
16, The calculation method according to claim 15, 
10 characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j,B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if primary mode, secondary mode and permutation 
15 signals are 0: 

si = el + e2 + e3 + e4 
s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 
20 2) if the primary mode signal is 0 and the 

secondary mode signal is 1: 
si = el + e4 
s2 = e2 
s3 = e3 
25 s4 = el - e4 

3) if the primary mode signal is 0 and the 
permutation signal is 1: 

si = (e3 + e4) - (el + e2) 
s2 = el - e2 
30 s3 = e4 - e3 

s4 = el + e2 + e3 + e4 
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4) if the primary mode signal is 1 and the 
permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
5 s3 = el - A.e3 + B.e4 

s4 = -e2 + B,e3 + A.e4 

5) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = el - A.e3 + B-e4 
10 s2 = -e2 + B.e3 + A.e4 

s3 = el + A,e3 - B.e4 
s4 = e2 + B,e3 A,e4 

17. The calculation method according to claim 10, 
characterized in that, in transformation step Ep, each 

15 calculation block comprises: 

2^-1 internal butterflies and a peripheral 
butterfly for the even values of index p, and 
2^-1 internal butterflies, otherwise. 

18. The calculation method according to claim 17, 
20 characterized in that the secondary mode signal is 1 if 

the peripheral butterfly is used for the first 
transformation step with p being even, and otherwise 0, 

19. The calculation method according to claim 18, 
characterized in that, for four input samples el, e2, 

25 eS, and e4, and for a complex coefficient W^=A+j .B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0 and the 
secondary mode signal is 1: 
30 si = el + e2 

s2 = el - e2 



e e E^B :L E! 7' IE 



, 0 i 7' Ol Z 



44 



s3 = e4 - e3 
s4 = e3 + e4 

2) if primary mode/ secondary mode and permutation 
signals are 0: 

5 si = el + e2 + e3 + e4 

s2 = el - e2 
s3 = e4 - e3 

s4 = (el + e2) - (e3 + e4) 

3) if the primary mode and secondary mode signals 
10 are 0 and the permutation signal is 1: 

si = (e3 + e4) - (el + e2) 
s2 = el - e2 
s3 = e4 - e3 
s4 = el + e2 + e3 + e4 
15 4) if the primary mode signal is 1 and the 

permutation signal is 0: 

si = el + A.e3 - B.e4 
s2 = e2 + B.e3 + A.e4 
s3 = el - A.e3 + B,e4 
20 s4 = -e2 + B.e3 + A,e4 

5) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A.e4 
25 s3 = el + A.e3 - B.e4 

s4 = e2 + B.e3 + A.e4 
20. The calculation method according to claim 8, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 
30 have 
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four inputs for receiving input samples and 
four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
5 inputs, 

in order to selectively apply different 
transformation operations to the input samples, each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 

10 and a coefficient admitted on said corresponding 
additional inputs , 

and in that the final step furthermore performs an 
addition and subtraction between the first and the last 
output sample provided in the last transformation step. 

15 21, The calculation method according to claim 20, 

characterized in that, in transformation step Ep, each 
calculation block comprises one peripheral butterfly 
and 2^-1 internal butterflies* 

22, The calculation method according to claim 21, 

20 characterized in that, for four input samples el, e2, 
e3, and e4, and for a complex coefficient W^=A+j ,B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode signal is 0: 
25 si = el + e2 

s2 = el - e2 
s3 = e4 - e3 
s4 = e3 + e4 

2) if the primary mode signal is 1 and the 
30 permutation signal is 0: 
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s2 = e2 + B.e3 + A.e4 
s3 = el - A,e3 + B.e4 
s4 = -e2 + B,e3 + A.e4 
3) if the primary signal is 1 and the permutation 
5 signal is 1 : 

si = el - A.e3 + B.e4 
s2 = -e2 + B.e3 + A,e4 
s3 = el + A.e3 -B,e4 
s4 = e2 + B.e3 + A.e4 
10 23. The calculation method according to claim 9 or 

characterized in that the first and second binary 
addresses of \i bits are generated for each butterfly, 
each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
15 address being greater than the first binary address. 

24. The calculation method according to claim 23, 
characterized in that said first and second binary 
addresses are consecutive and an internal butterfly is 
involved. 

20 25. The calculation method according to claim 23 

or"'''^^, characterized in that, if a peripheral butterfly 
is involved, the p+2 low-order bits of the first 
address are equal to 0, and the p+2 low-order bits of 
the second address form a number equal to 2^'^"''-l, 

25 26. The calculation method according to claim 24 

or ^"^^ characterized in that the address of the two 
other samples to be applied to the inputs of the 
butterfly, be they peripheral or internal, are obtained 
by inverting the (p+2) low-order bits of said first and 

30 second produced addresses. 
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27. The calculation method according to claim 26, 
characterized in that even-numbered address samples and 
odd-numbered address samples are stored in two separate 
memories . 

5 28. The calculation method according to claim 25, 

characterized in that the value of the parameter s of 
the coefficient assigned to an internal butterfly in 
transformation step Ep is coded by p-2 bits, and is: 

if p+l=vi-2, the number formed by the p+1 
10 low-order bits of the second binary address produced 
for said internal butterfly, 

if p+l<vi-2, the number formed by the p+1 
low-order bits of the second binary address produced 
for said internal butterfly, followed by p-p-3 zero 
15 bits at the end of the number, 

if p+l>p-2, the number formed by the p+1 low- 
order bits of the second binary address produced for 
said internal butterfly, minus its p-p-1 low-order 
bits . 

20 29, The calculation method according to claim 4, 

^in turn dependent on claim 3, in turn dependent on 
claim 2,^ characterized in that in each transformation 
step Ep, the butterflies are distributed among 2^ 
calculation blocks, 

25 in that each calculation block comprises one 

peripheral butterfly and N/2P''^-1 internal butterflies, 

in that the peripheral butterfly of the rank a 
calculation block in transformation step Ep transforms 
the input samples of rank 2^-Pa, 2^-*^a+2^-'*-^-l , 2^-W2^-P- 

30 ^, 2^"^a+2^"^-l into output samples of the same rank, 




and in that the internal rank x butterfly of the 
rank a calculation block in transformation step Ep 
transforms the input samples of rank 2^'^a+2x+l/ 2^" 
Pa+2x+2, 2^-Pa+2^-P-^-2x-3, 2^-*^a+2^-P-2x-2 into output 
5 samples of the same rank . 

30. The calculation method according to claim 29, 
characterized in further comprising a final step of 
modifying the sequence of the output samples provided 
in the last transformation step and classifying them in 

10 ascending order of index n. 

31. The calculation method according to claim 29 
oiT^^S^ characterized in that each butterfly is assigned 
a coefficient W^, whereon the calculation inside the 
butterfly is based, said coefficient being equal to e" 

-1^5 j{2«s/N) ^ith s G [0.. N/4-1] for a fast Fourier transform 
and is equal to e^ with s g [0.,N/4-1] for an 

inverse fast Fourier transform. 

32. Calculation method according to claim 31, 
characterized in that the internal rank x butterfly of 

20 the rank a calculation block in transformation step Ep 
is assigned coefficient with 5 = (x+1) .2*^. 

33. The calculation method according to claim 32, 
characterized in that the butterflies for implementing 
the transformation steps are all of the same type and 

2 5 have 

four inputs for receiving input samples and 
four outputs for providing output samples, 

four additional inputs, respectively primary 
mode, secondary mode, permutation, and coefficient 
30 inputs. 
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in order to selectively apply different 
transformation operations to the input _^ samples , each 
operation being determined by the values assigned to 
the primary mode, secondary mode, permutation signals, 
5 and a coefficient admitted on said corresponding 
additional inputs . 

34. The calculation method according to claim 33, 
characterized in that, for each butterfly, the primary 
mode signal is 0 for a peripheral butterfly and 1 for 

10 an internal butterfly, 

in that the permutation signal is 0 for the even 
rank calculation blocks, including rank 0, and 1 for 
the odd values. 

35. The calculation method according to claim 31 
15 or'^^^ characterized in that the secondary mode signal 

is 1 if the butterfly, be it peripheral or internal, is 
used for the first transformation step, and otherwise 
0. 

36. The calculation method according to claim 35, 
20 characterized in that, for four input samples el, e2, 

e3, and e4, and for a complex coefficient W^=A+j.B, the 
butterfly delivers the following output samples si, s2, 
s3, and s4 

1) if the primary mode and secondary mode signals 
25 are 0: 

si = (el + e2)/2 
s2 = (el - e2) /2 
s3 = (e4 - e3) /2 
s4 = (e3 + e4) /2 
30 2) if the primary mode signal is 0 and the 

secondary mode signal is 1: 
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si = [ (el + e4) /2-e2] /2 
s2 - [ (el+e4) /2-e2] /2 
s3 = [e3- (el-e4) /2] /2 
s4 = [e3+ (el+e4) /2] /2 
5 3) if the primary mode signal is 1 and the 

permutation signal is 0: 
si = (el+e3)/2 
s2 = (e2+e4) /2 

s3 - [(el-e3).A - (e2+e4).B]/2 
10 s4 = [-(el-e3).B + (e2+e4).A]/2 

4) if the primary mode signal is 1 and the 
permutation signal is 1: 

si = [ (el-e3) .A - {e2 + e4) .B] /2 
s2 = [-(el-e3).B + (e2+e4).A]/2 
15 s3 = (el+e3) /2 

s4 = (e2-e4)/2 

37. The calculation method according to claim 33, 
characterized in that the first and second binary 
addresses of ]i bits are generated for each butterfly, 

20 each binary address corresponding to the rank of an 
input sample of said butterfly and the second binary 
address being greater than the first binary address. 

38. The calculation method according to claim 37, 
characterized in that said first and second binary 

25 addresses are consecutive and an internal butterfly is 
involved. 

39. The calculation method according to claim 37 
or characterized in that, if a peripheral butterfly 
is involved, the ]a-p low-order bits of the first 

30 address are equal to 0, and the u-p low-order bits of 
the second address form a number equal to N/2^'*"^-l. 
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40, The calculation method according to claim 38 
or """S^ characterized in that the address of the two 
other samples to be applied to the inputs of the 
butterfly are obtained by inverting the p-p low-order 

5 bits of both produced addresses. 

41. The calculation method according to claim 40, 
characterized in that even-numbered address samples and 
odd-numbered address samples are stored in two separate 
memories . 

10 42. The calculation method according to claim 41, 

characterized in that the value of the parameter s of 
the coefficient assigned to an internal butterfly in 
transformation step Ep is coded by vi-2 bits, and is: 

if p-p-l=u-2, the number formed by the p-p-1 
15 low-order bits of the second address produced for said 
internal butterfly, 

if ]i-p-l<]jL-2, the number formed by the p-p-1 
low-order bits of the second address produced for said 
internal butterfly, followed by p-1 zero bits at the 
20 end of the number, 

if vi-p-l>p-2, the number formed by the u-p-1 
low-order bits of the second address produced for said 
internal butterfly, minus its p+1 low-order bits. 
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